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Abstract'^ 

Recently we showed that a simple model of network rewiring could be solved exactly for any time and 
any parameter value. We also showed that this model can be recast in terms of several well known models 
of statistical physics such as Urn model and the Voter model. We also noted that it has been applied to 
a wide range of problems. Here we consider various generalisations of this model and include some new 
exact results. 



1 Introduction 

Graphs with a constant number of edges and vertices but which evolve by rewiring those edges are a classic 
network model as exemplified by Watts and Stogatz [1] ([2, 3] provide further examples). Such network 
evolution may also be recast as other types of statistical physics models (for example see [4, 5, 6, 7]). 
As many real systems are effectively of constant size, non growing networks can also be used to model 
a wide range of data: the transmission of cultural artifacts such as pottery designs, dog breed and baby 
name popularity (as in [8, 9, 10, 11]), the distribution of family names in constant populations, and the 
diversity of genes. In this paper we look at various extensions to a model of network rewiring for which 
an exact solution [12] was presented at ECCS06 [13] and in more detail in [14]. 



2 The Model 

We will study the rewiring of a bipartite graph consisting of E 'individual' vertices connected by one 
edge only to any one of N 'artifact' vertices, as shown in Fig. 1. At each time step two choices are made. 
With probability Uj^ an individual is chosen. It is the artifact end of its single edge, connected to the 
'departure' artifact, which is to be rewired. An 'arrival' artifact is also selected with probability Ua- 
Only after the choices are made is the network altered by rewiring the chosen edge so that its artifact 
end is moved from the departure to the arrival artifact. Note we do not explicitly exclude the possibility 
that the departure and arrival artifacts are the same. The individual vertices always retain one edge 
while the degree k of the artifact vertices is changing in time, only its average degree (fc) — E/N is 
constant. It is the distribution of the artifact vertices at time n{k,t), and its probability distribution 
p{k,t) = n{k, t)/N , that we study. This process can be viewed in many other ways [14]. 

The evolution of the degree distribution in the mean field approximation is described by the master 
equation [12, 13, 14] 
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Fig. 1. The bipartite graph has E 'individual' vertices, each with one edge. The other end of the edge is connected 
to one of A'^ 'artifact' vertices. If the degree of an artifact vertex is k then this artifact has been 'chosen' by k 
distinct individuals. At each time step a single rewiring of the artifact end of one edge occurs. An individual is 
chosen (number (2i — 1) here) with probability 11 r which gives us the departure artifact (here D). At the same 
time the arrival artifact is chosen with probability 11 a (here labelled B). After both choices have been made the 
rewiring is performed (here individual (2i — 1) switches its edge from artifact D to B). 

n{k,t+ 1) -n{k,t) 

= n{k + 1, t)i7fi(fc + 1, t) (1 - nA{k + l,t)) 

-n{k, t)nR{k, t) (1 - HAik, t)) - n{k, t)nA{k, t) (1 - Hnik, t)) 

+n{k ^ l,t)nAik - l,t) {I - nR{k ^ l,t)) . (1) 

For our physical problem the removal probability must always satisfy Ufiik = 0) = and Uiiik ^ E) = \. 
In addition for physical solutions we must have n{k,t) = if k<0 or k>E. The presence of the 
factors of (1 — 77) ensure that if the degree distribution initially satisfies its physical boundary condition, 
n{k, t — 0) — ii k < 01 k > E, then this boundary condition is automatically satisfied at all times'*. 
The factors of (1 — 77) are not seen in the master equations of the literature [2, 3, 4, 15] and correspond 
to events where the arrival and departure artifacts are chosen to be the same^. 

In general the master equation (1) gives the evolution only in the mean- field approximation because we 
are taking an ensemble average over many instances of the stochastic evolution and using the product of 
averages where we should have the average of products. However when the normalisations of probabilities 
77^ and 77^ are constant then the master equation (1) may be exact. The most general 11 and Ua for 
which this is true is 

"r^^, nA^Pr^ + a-Pr)^, iE>k>0). (2) 

We will restrict ourselves to these forms and therefore our analytic results are exact [14]. Thus we are 
choosing our arrival edge with a mixture of preferential attachment (probability (1 — Pr)) and random 
attachment (probability Pr). The removal artifact is found by choosing the artifact end of a randomly 
selected edge — 'preferential removal'. The use of probabilities proportional to k can emerge naturally 
through short range searches of many networks, since the probability of arriving at a vertex on a random 
graph is proportional to its degree [16, 17, 18]. 

Not only is the master equation exact for our chosen probabilities (2) but the exact solution for the 
degree distribution n{k,t) may be found for any finite parameter value. This may be done in terms of 
(E + 1) eigenfunctions a;'^™^(fc) and their corresponding generating functions G'^™^(a;) 



* For this to be true it is absolutely vital that we have the factors of (1 — i7ij(fc)) to ensure that with the 

condition Tliiik = 7?) = 1 we do not include processes where an artifact with E edges is lost because we are 

adding another edge (the third term in (1) for k = E). 
^ These events occur with probability {null a)- Since the network is unchanged by such events, we must exclude 

such events from the evolution of n{k,t) and the factors of (1 — 77) implement this. It is an approximation to 

drop these terms. In our model this is not justified for certain parameter values. 
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E 

n{k,t)^Np{k,t) - ^c™(A™)*c^('")(/c), (3) 



m=0 
E 



G{x,t) := ^a;'=n(fc,t) = ^ c„(A,„)*G(")(x) , (4) 

k=0 TO=0 

E 

G('")(a;) ^a;'=tj('")(fc). (5) 

fc=0 

The solution is found to consist of simple combinations of the Hypergeometric function F{a, b; c; x) [14] 
G(™^(x) = (1 -a;)'"F(a + m,6 + TO;c;a;) (6) 

„n .„,yr r(a + m + or(6 + m + or(c) , 

^ r(a + m)r(6 + m)r(c + Z)(;!) 

with corresponding eigenvalues, 

A„ ^ 1 - - m(m - 1) , S>m>0. (8) 

The eigenvalues satisfy A,„ > A^+i except for = when Aq = Ai = 1. 

There is a unique long time equilibrium distribution which can be of one of two phases, as Fig. 2 
shows. For pr < we get a condensate, most individuals attach to a single artifact. For pr S> E^^ 
we get a power law degree distribution of unit slope, with an exponential cutoff n(k) « k~^e~^^ [C^ = 
— In (1 — Pr))- There is a smooth transition between the two except in the i? — > oo limit. 




Fig. 2. Plots of the degree probability distribution 
function p{k) = n{k)/N ior N = E = 100 and var- 
ious Pr = 1 (red crosses), 10/_E (green circles), 1/E 
(blue stars) and 0.1/-E (magenta squares). Note that 
Pr = 1/-E is almost a pure power law for all values of 
k. 



Using the fact that the number of artifacts TV and the number of edges E are constant gives cq = N 
and Ci = so that the cigenmode numbered one (m = 1) never contributes. Thus the approach to 
equilibrium of most quantities occurs on a timescale 

r2 = -[ln(A2)]-i (9) 

as illustrated in Fig. 3. This means that if we have rewired most of the edges once and almost never 
used random attachment, i.e. Pr ^ E~^, then the approach to equilibrium is slow, T2 = 0{E^). However 
for other cases, Pr ^ E~^ , the small amount of randomness gives a rapid approach to equilibrium after 
every edge has been rewired just a few times. The initial conditions determine the remaining Cm (fn > 1)- 

Of particular interest are the homogeneity measures Fn (t) which are the probability that n randomly 
chosen but distinct edges all share the same artifact. These are given by 
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Fig. 3. Plots of p{k) from simulations (data points) and the exact analytic results (lines) for E — N = 100 
with pr = 10/i? on the left and p,. = 0.1/E on the right. The results are shown at four different times: t k, T2 
(red, crosses), t ~ 2t2 (green, circles), t « 3r2 (blue, stars) and to equilibrium (magenta, squares). The initial 
configuration has one edge per artifact. The data points are averages over 10^ runs while the lines are the exact 
analytic results. 
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(10) 



The properties of the Hypergeometric function mean we can express the solutions in terms of fixed 
fractions of a large number of fixed Gamma functions with all the dependence on time and the initial 
conditions is carried by factors of Cm(Am)*- Also the n-th homogeneity measure Fn{t) has contributions 
only from the eigenfunctions m < n. 

For instance the most useful homogeneity measure is F2{t): 



F2{t) = ^^2(00) + (A2)* (F2(0) - ^^2(00)) , ^^2(00) 



l+p,.((fc) -1) 

l+Pr{E-l) ' 



(11) 



while the initial conditions set ^2(0). 



3 Phase Transitions of Unipartite Graphs in Real Time 

The construction of MoUoy and Reed [19] gives a unipartite graph of a given degree distribution but is 
otherwise random. In terms of our bipartite graph this is equivalent to taking pairs of individual vertices 
and merging the edge ends ('stubs') coming out of these individual vertices. The individual vertices are 
then thrown away. This is illustrated in Fig. 4. 

The rewiring of our bipartite model is then equivalent to a rewiring of the projected unipartite graph 
with the same linear attachment and removal probabilities, also illustrated in Fig. 4. Since the degree 




Fig. 4. One projection of the bipartite graph of Fig. 1 onto this undi- 
rected unipartite graph. Let a{i) be the artifact vertex connected to the 
individual vertex i in the bipartite graph. Then we take pairs of individual 
vertices (2i) and (2i — 1) in the bipartite graph and connect their associ- 
ated artifacts a{2i — 1) and a[2i) in the undirected graph. The rewiring 



event of Fig. 1 now become a rewiring of the (D,H) edge to a (B,H) edge. 



distribution of our artifact vertices is also the degree distribution of the unipartite graph, all our results 
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can be applied directly to such graphs. For instance for pr = I we capture the degree distribution of the 
original Watts and Stogatz model^ [1]. 

Analytic expressions for the global properties of such random graphs in the infinite A'^ limit depend 
on a ratio, z, of the second and first moments of the degree distribution [16, 17, 18, 19] 

zit):=^-^-l = {E-l)F2it). (12) 

There is a phase transition in the properties of such infinite random graphs at z = 1. This occurs when 
there is one tadpole (an edge connected at both ends to the same vertex) in the unipartite graph. In 
particular for z > 1 the average distance between two vertices in the giant component, {I), may be 
estimated to be'' [18] 

ffl = =«a±M£H7s + |, ,,.o.„r2. (13, 

ln(2;) 2 

For simplicity we consider graphs where N = E which start with each artifact connected to only 
one individual so n{k,t = 0) = E5k,i- The projected unipartite graph has (k) = 1 and initially ^2(0) = 
2;(0) = 0. li Pr ^ 0{E~^) then the equilibrium configuration is reached quickly in t ~ 0(t2) = 0{E) 
steps. Only when we start to get a high degree node, so a condensate is forming and Pr < 0{E~^), do 
we get a slower approach to equilibrium on a time scale T2 = 0{E'^). The phase transition in infinite 
random graphs occurs at z = 1. In our case, our projected graphs start from z(0) = but they reach 
z = 1 very quickly at ti « E/2 unless (1 —Pr) '> 0(E^^). That is even if the evolution to the equilibrium 
distribution is slow, provided a reasonably large degree node exists, i.e. there is significant amount of 
copying, a large component emerges in the projected unipartite graph quickly, typically at ti » E/2, 
since 

ln(l -((£;- 1)^2(00))-^) .... 

*i = 1 N , (14) 

ln(A2) 

E E , , 1 

(l-p.)»-. (15) 



2{l+pr{{k)-l)yE-l' ' ^•^'"-^ E- 

The numerical results for the evolution of the properties of the projected unipartite graph are shown 
in Fig. 5. The parameter z reaches the value 1 at (ti/E) « 0.5 ± 0.0002 as expected. This is close to, 
but not exactly equal to, (tp/E) — 0.535 ± 0.005, where tp is the time at which the average distance and 
diameter of the largest c;omponent peak. The second derivative in time of the number of vertices in the 
largest component also suddenly switches sign at exactly the same time tp. The value of z at this time 
is z{tp) = 1.06 ±0.01. 

Motivated by the approximate expression for the distance in the largest component of a large random 
graph (13), we find that the inverse distance for the parameters used in Fig. 5 is well fitted by the form 
aln(z — 0.06) + b + cz + dz^ but with different values either side of the peak time^. The fit is shown in 
Fig. 6. The deviations from the predicted z = 1 transition point seem to be finite size effects. The peaks 



Strictly speaking we choose a random edge to rewire while Watts and Stogatz [1] rewired in a systematic 
manner. 

This formula must be adapted from [18] to take account of the existence of vertices of zero degree. Analytic 
derivations of such global properties use an ensemble of graphs over which there is always a finite probability 

of getting from any one vertex of degree fci > to a vertex of degree kj > in a finite number of steps. In any 
one graph this need not be true. On the other hand numerically we measure the average distance in the largest 
component of one graph considering only vertices in the largest component. We then average this result over 
the ensemble of graphs. Numerical evidence suggests that this numerical measurement has the same qualitative 
behaviour as the analytic formula. 

For early times, z < z(tp), we have a = -0.107 ± 0.006, b = 0.30 ± 0.02, c = -0.42 ± 0.04 and d = 0.14 ± 0.03 
(fit excluded the four points with lowest z values) while for late times and z > z{tp), we have a = +0.85 ±0.02, 
b = 0.019 ± 0.02, c = -0.002 ± 0.06 and rf = 0.0008 ± 0.0003. These fits have R'^ = 0.9995 and R'^ = 0.9999 
respectively and errors are at 95% confidence level. However a polynomial works almost as well, at least near 
t — tj) . 
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Fig. 5. Properties of the undirected random graph 
formed using a Molloy-Reed type projection [19] . The 
underlying bipartite graph has N = E = 10^ start- 
ing from F2(0) = and rewired with pure copying 
{pr = 0.0). Results are calculated for each instance 
and then averaged over a total of 1000 runs. 
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Fig. 6. The inverse distance of the largest compo- 
... ■ " nent for the same projected networks as the previous 

if figure. The points are the data (errors are smaller 

than the symbol size) and the lines are the best fits 
< < < < < < J to the form aln(z — 0.06) + b + cz + dz^ . The lower 

0.5 1 1.5 2 2.5 3 3.5 4 ^ ' 

z figure shows the residuals illustrating the good fit. 



« -2- 



are sharper and closer to z = 1 as the network get larger^ but with the same (/c), pr and ^2(0) — 0. 

The network shown is evolved with pure copying pr = Q so in this case the equilibrium distribution, 
a complete condensate F2 = 1, will emerge only on a long time scale of T2 = — ln(l — 2E~'^) ~ 0{E^). 

One way to look at this transition is to use the interpretation of the model in terms of cultural 
transmission [8, 9, 10, 11, 12, 13, 14]. In this case the bipartite graph represents individuals who are 
choosing artifacts by either copying the choices made by another individual (preferential attachment) or 
by making their own innovation (random attachment). Suppose we now imagine that each person has 
two copies of an artifact. The unipartite graph is then one expression of the relationship between objects 
as defined by the choices made by individuals. For instance one could imagine asking people to categorise 
their two favourite pairs of shoes and each artifact could represent a different category, e.g. one artifact 
might represent black leather lace up shoes. The unipartite projection gives a metric in artifact space 
as defined by the choices made by the individuals. The phase transition in the unipartite network then 
marks the point where the individuals have reached some sort of consensus as the artifacts now form a 
Giant Connected Component given the metric provided by the individuals' choices. 



4 Voter Models and Individual Networks 

One possible generalisation of our rewiring model is to add a second graph connecting the individual 
vertices which we will call the Individual graph. When an individual rewires using preferential attachment 

^ FoT N ^ E and Pr = 0.0, we find: iV = 10^ tp/E = 0.66 ± 0.04; N = 10*, tp/E = 0.57 ± 0.01; iV = 10^ 
tp/E = 0.535 ± 0.005. Estimated from an ensemble of 1000 independent runs for each value oi N = E. 
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they copy the artifact chosen by one of their neighbours in the Individual network. With pr = and 
= 2 we obtain the basic Voter model [6, 7]. Our model corresponds to having a complete graph for 
the individual network but with the addition of a random rewiring process, Pr > 0, and an arbitrarily 
large number of choices, N > 2. Neither of these cases is studied in the Voter model literature where 
the focus is on different types of individual networks and any analytic results are only available for the 
E ^oo limit [7]^° 

H rfi „ 



Fig. 7. Equilibrium artifact degree distribution p{k) for dif- 
ferent Individual graphs of 100 vertices and average degree 4: 
Erdos-Reyni (red pluses), Exponential (green circles), Barabasi- 
Albert (purple squares), periodic lattices of two (grey crosses) 
and one (blue diamonds) dimension. The line is the analytic re- 
sult for a complete Individual graph while the other results are 
taken over an ensemble of 10^ Individual graphs. N = E = 100, 

Pr = 1/E. 

Results for the equilibrium distribution show that it is qualitatively unchanged by the type of indi- 
vidual graph^^ except for the case of a one dimensional ring, as shown in Fig. 7. 

We will use two quantities to study the behaviour of the model. Our quantity F2 of (10) is a measure 
of the global homogeneity. An equivalent measure which takes account of the local properties of the 
Individual network is the average interface density, (p), the probability that any two individual vertices 
connected by the Individual graph have a different artifact. If the graph is complete or if the Individual 
graph is ignored (pr — 1) then from (11) we have {p)t = I — F2{t). Otherwise for two reasons we expect 
that {p)t ^ 1 — F2{t) and that both would both differ from value obtained for a complete Individual graph 
as derived from (11). First because of the explicit reference to the Individual graph in the definition of 
{p)t but not in ¥2- Second, the structure imposed by the Individual graph will, in general, effect both 
the evolution timescale and, for pr > 0, equilibrium degree distributions as compared to the complete 
Individual graph case. 

We can see these differences if we compare the equilibrium values reached on lattices of different 
dimensions but with some randomness present (otherwise {p)t — (1 — P2{t)) because both are zero). As 
Fig. 8 and tabic 1 show, the local and global homogeneity measures {p) and (1 — F2) are close to the 
analytic result for large dimension lattices with short network distances. As we take lattices of smaller 
dimension, F2 gets much larger than the analytic result, and {p) much smaller. Table 2 shows a similar 
effect as we increase pr- 

The time scale of the approach to consensus is often studied in Voter models. If pr > 0, so there 
is no absolute consensus, the approach to equilibrium, as measured by F2 and p, is controlled solely 
by the time scale of the second eigenvalue T2 of (9) if the Individuals are connected by a complete 
network. For general Individual networks the evolution of the global F2 or local (p) takes the same form 
as (11), aexp(— t/^o) + c. However, as one might expect, the formation of small patches of consensus 
between nearest neighbours, measured by (p), happens faster than the emergence of a global consensus. 

For Pr > we can think of our model as including two graphs. The first, as mentioned above, is a graph 
connecting Individuals. Preferential rewiring is done by performing a random walk of length one on this graph, 
and copying the choice of the resultant Individual. The second graph connects the Artifacts [14] . If this graph 
is a complete graph (with tadpoles) then a random walk on the graph gives the random attachment pr term 
appearing in TTa of (2). One may imagine many practical problems where the Artifact network is not so trivial. 
For this paper, however, we only consider the case of a complete Artifact graph. 

In this article, the lattices are cubic (Z*^). In Fig. 7 next-to- nearest and nearest neighbours are connected in the 
one- dimensional ring. In all other cases the only nearest neighbours connected in the lattice Individual graphs. 
The Exponential and Barabasi-Albert graphs are connected graphs Individual graph degree distributions of 
Pind(fc) oc exp{— ^fc} and Pind(fc) cx: [fc(fc + l)(fc + 2)]^^ respectively. The results also support the claim in 
[12, 13, 14] that some results for a Minority game played on an Erdos-Reyni graph can be understood in terms 
of our results for the rewiring model. 
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Fig. 8. Homogeneity measures for various lattices 
against t/E. The black solid line represents the an- 
alytic 1 - F2{t) for N = 2, pr = l/E a.nA E = 729. 
Numerical results for 1 — F2{i) (triangle highlights) 
are plotted for 1-d (red), 2-d (purple) and 3-d (blue) 
regular lattices. The average interface densities (p) 
are plotted as circles. Averaged over 1000 runs. 



Dim 


to{F2)/T2 


to{p)/T2 


1-F2(oo) 


p{oo) 


Id 
2d 
3d 


1.25 (3) 
1.19 (1) 
1.10 (1) 


0.0241 (4) 
0.74 (1) 
1.06 (1) 


0.47466 (2) 
0.3494 (1) 
0.2898 (1) 


0.0261 (1) 
0.1558 (1) 
0.2120 (1) 



Table 1. Table of time scales and the limiting value of the evolution of 1 — F2 and p where the Individual graphs 
are periodic lattices with nearest neighbour connections only. The complete graph has 1 — F2{t = 00) ~ 0.25017 
and T2 ~ 1.3295*. Extracted from the data of Fig. 8 by fitting to aexp(— t/io) +c with the estimated error in the 
last digit give by the numbers in brackets. 

as measured by F2 ■ This is accentuated if there is a large distance between individuals as the comparison 
between lattices of different dimensions in Fig. 8 and in table 1 show. Varying also shows that local 
equilibration is faster than global but there does seem to be a marked difference between <C 1/E 
and Pr > 1/E as table 2 shows. For prE <^ 1, local equilibration is a little slower than occurs on 
complete graph. However for prE > 1, this randomness brings local equilibrium an order of magnitude 
faster than was the case with a complete graph. It shows that a little bit of randomness can speed up 
local equilibration but not if an overwhelming consensus is going to emerge. 



Pr/E 


T2 


to{F2)/T2 


to{p)/T2 


1 — ^2(00) exact 


1-F2(CX)) 


p{oo) 





79999 


2.131 (6) 


1.8 (1) 





0.0003 (3) 


0.003 (3) 


0.1 


72743 


1.988 (5) 


1.81 (2) 


0.04546 


0.0827 (3) 


0.0403 (4) 


1 


40050 


1.34 (3) 


0.15 (2) 


0.25031 


0.3375 (3) 


0.167 (1) 


10 


7289 


1.01 (6) 


0.145 (3) 


0.45558 


0.47433 (4) 


0.2634 (1) 


100 


794 


0.7 (2) 


0.40 (1) 


0.49628 


0.49712 (1) 


0.3813 (1) 



Table 2. Table of time scales in units of T2 and limiting value in units of the exact value for c = F2{t = 00) for 
400 individuals connected by a square lattice. Data averaged over 5000 runs for pr = and Pr = 1/E, and 1000 
runs for all others. It was fitted to aexp{—t/to) + c. 



In table 3 we see that the time scales for the exponential decay obtained from fitting our data are 
roughly in line for the predictions made for the completion time in this model [6, 7, 20] on a lattice, 
to ~ 0{E) for a one-dimensional lattice, to ^ 0(ln(£')) in two dimensions. However some discrepancies 
suggest more work is needed. 

The main result to draw from table 3 is that the evolution towards equilibrium, its time scale and 
final value, are independent of the number of artifacts. This is to be expected at small Pr given our 
linear attachment probabilities as this gives our model certain scaling properties [14]. Suppose we have 
the consensus emerging picking out one of our N artifacts and we merge the remaining — 1 artifacts 
into one artifact. The probability of an individual copying the consensus artifact or one of the remaining 
artifacts is exactly the same as if we had a model with N ^ 2 and the same (1 —pr)- The only difference 




See also Fig. 10. 
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N = 


2 


= 10 


Dim 


E T2 


to{l - F2)/t2 


to{p)/T2 


to(l - F2)/t2 


to{p)/r2 




100 4999.5 


18.6 (1) 


0.42 (2) 


18.9 (1) 


0.41 (1) 


Id 


200 19999 


37.9 


0.21 (1) 


39.3 


0.21 (1) 




400 79999 


76.3 


0.18 (1) 


75.1 


0.18 (1) 




1000 500000 


115.7 


0.069 (2) 


137.8 


0.070 (3) 




400 79999 


2.131 (3) 


1.8 (1) 


2.109 (2) 


1.8 (1) 


2d 


900 405000 


2.3 


2.1 


2.3 


2.1 




2500 3130000 


2.8 


2.1 


2.8 


2.1 



Table 3. Table of time scales in units of T2 for E individuals connected by a one- or two-dimensional torus 
found by fitting the data to aexp(— f/fo) + c. For two and ten types of artifact, pr = 0. Data was averaged over 
1000 runs for the largest lattices down to 50 runs for the smallest lattices. Where the error is known reliably, the 
numbers in brackets specify the error in the last digit. 



is that when a random innovation event occurs, with probability p^, the non-consensus artifacts are 
preferred to the single consensus artifact by a factor of [N — 1). Thus for TV ^ 2 the random events 
are more likely to destroy the emerging consensus than in the Voter model but only if Pr ^ 0. For 
the extreme case of = we see the expected lack of dependence on N in table 3. The only effect of 
increasing the number of artifacts in our results comes from starting from a homogeneous initial condition 
so F2 = which is further away from F2 — 1 and consensus, see Fig. 9. 
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Fig. 9. Plots of the analytic (1 — F2{t)) (hues) and (p)t for 
1-d periodic lattices against t/E for pr = G and E = 100 
(red, far left), 200 (purple), 400 (blue) and 1000 (black, far 
right). The lower circles and the solid lines represent N — 2 
while the higher squares and dashed lines are for A'^ = 10. 
Data are averages over 1000, 1000, 500 and 100 runs for 
increasing lattice sizes respectively. 



5 Two Types of Individual 

Another variation of our original model is to introduce two types of individual, labelled X and Y . At 
each time step we first pick which type of individual to update; with probability we select at random 
one the X-type individuals. We rewire its artifact end, choosing its arrival artifact in one of three ways: 
at random, by copying the existing choice of one its own type of individual, or finally copying the existing 
choice made by a random individual of the opposite type. These arrival probabilities may be different 
for the two types so we have four independent arrival probabilities and one departure probability. Add 
in the freedom to choose different numbers of X and Y individuals, and Ey, and N the number of 
artifacts, we find we have eight free parameters. The degree distribution is now n(kx, ky;t), the number 
of artifact vertices which have (ky) edges to X (Y) type vertices at time t. 

The question is can we solve this system analytically? By keeping our probabilities linear in degree 
and because our normalisations are constants of the evolution, our mean field equation is again exact, for 
the same reasons as in the original model [14]. Writing in terms of the generating function G{x,y,t) :— 
Sfj=o Sf'J'=o x''''y''^n{kx,ky, t) we find that we can again spfit this into {E^ + l){Ey + 1) eigenfunctions 
which we label with a pair of indices (M, A): 
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■.= f^Y.i^-iy(y-iyfl^^\ (16) 

where f^^^^ are constants. The eigenfunctions satisfy a two-dimensional second order PDE. We have 
not found a full solution but we can reduce this to a one-dimcnsional problem. Since we express our 
eigenfunctions in powers of (a; — 1) and {y — 1), the eigenfunctions satisfy^'^ fij ^^^ = if i + j < M 

for any integer < M < + Ey. At the same time the equations for the coefficients f-^^^ and 
eigenvalues Xma, where i + j = M, involves no coefficients where i + j > M. Thus the label A indexes 
the allowed values of i and j given the constraint i+j = M. Finding the eigenvalues is therefore a matter 
of solving a set of {inm{Ex,Ey) + 1) linear equations. This also gives the coefficients of the eigenfunctions 
for i + j = M. The remaining values i+j > M may be found iteratively though this is generally a 
non-trivial problem. 

However, we have seen that much information is encoded by the first and second moments of the 
degree distribution. The general homogeneity measures are given by 



_ r{E^ + l) r{Ey + l) d^+^G{x,y;t) 



r{E^ -m + l) r{Ey - n + 1) aa;™5y" 



(17) 



So wc only need the three second order homogeneity measures, m + n = 2. These have contributions 
only from the i+ j <2 coefficients and therefore only the M < 2 eigenfunctions contribute. The system 
of equations for such coefiicients reduces to solving for the eigenvalues and eigenfunctions of a three 
dimensional system, which has an exact, if lengthy, algebraic solution. The basic results though arc that 
the only equilibrium solution is given by the single M = cigcnfunction where Aq = 1, f^Q = N, 
/}q^ — E^, /q"-' = Ey, and the three coefficients , f[^\ satisfy 



f(0) 



(f2o\ ( a,{l + a,){E, - I) 

f[f = T-M a^{d^ - a^/E^)Ey + ay{dy - ay/Ey)E^ I , (18) 

\/(°7 V ay{l + ay){Ey-l) 

-2a^(l + E-^) -a^d^il - E^) 

-2aydy axPx + ay(3y -2axdx | , (19) 

-aydy{l-Ey) -2ayil + Ey)^ 

" -(A:.) , d, = , a, = , /3x = Pp--\ (20) 



Ppxx Ppxx Ppxx 

Switching labels (a; y) gives the similar y subscript parameters. 

The M = I eigenfunctions again give no contribution to any physical quantity since the first moments 
are constant. The second moments are given in terms of the lowest coefficients of one of three M = 2 
eigenfunctions, the /2o'*\ (assuming Ex,Ey > 1), which satisfy 

(A2A - l)f('^^ = Tf(2^) , f(^^) = /(f ^) . (21) 

There is a large parameter space to investigate but there are a few obvious limits. First one can 
scale the probabilities in proportion to the number of edges of each type so Qa = Ea/E, pra = Pr, 
Ppab = (1 —Pr)Eb/E where a, 6 G {x, y}, E = Ex + Ey. One can sec then that the total degree distribution 
given by G{x, x) is exactly as we had in the single type model. However we can now investigate the 
'chemical equilibrium' as the distribution of X and Y types, given by derivatives of G{x, 1) and G{l,y) 
respectively, will evolve differently if the initial conditions are different for each type. Another simple 
example is where Pxx = Pyy = which encodes the "complete bipartite graph example" of [7]. Our 
method allows one to extract exact expressions for the whole time evolution, not just order of magnitude 
estimates for the equilibration time. 



We use this to define our label M. 
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6 Conclusions 

In this paper we have looked at a variety of extensions to the basic network rewiring model of [12, 13, 14]. 
Studying the projection onto a unipartite graph gives us exact expressions for the time evolution of a 
finite sized system through a transition. 

We have also shown that adding an Individual network leaves the qualitative behaviour of the model 
is unchanged in terms of F2 . However quantitative differences are highlighted by comparison against the 
case of a complete graph for which our previous analytic work [12, 13, 14] provides exact analytic results. 
What we learn from this model is that the consensus (the condensate) may not be perfect, 1 > prE > 0, 
and it may emerge very slowly T2 ~ 0{E'^), but an effective consensus is always reached very quickly 
ii ~ 0{E). 

Finally we have shown how some progress can be made on solving models with more than one type 
of edge. In particular we show how the various homogeneity measures Fmn niay be found exactly. 
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Additional Figures 



These were not included in the proceedings. Data in some of the tables was derived from these curves. 




Fig. 10. Plots of 1 - F2(t) (triangles) and (p) (cir- 
cles) on a 2-d periodic lattice with side L — 20 
and N — 2 with attachment probabilities pr — 
(red), Pr — 0.1/ E (purple), pr = 1/E (blue) and 
Pr = 100/i5 (black). Solid coloured lines are the 
equivalent analytic results for 1 — F2{t) on a com- 
plete graph. Averaged over 5000 runs for pr = and 
Pr ~ 1/E and 1000 runs for all others. Data used for 
table 2. 
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Fig. 11. Analytic 1—F2{t) (solid lines) for a complete 
graph compared against numeric (p) (circles) on a 2- 
d periodic lattice with N = 2 and pr = 0. L = 20 
(red), i = 30 (blue) and L — 50 (black). Averaged 
over 5000, 500 and 50 runs respectively. Data used 
for table 3. 
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Fig. 12. Plots of (p) on a 2-d periodic lattice with 
Pr = for N — 2 (circles) and A'' = 10 (squares). 
L = 20 (red), L = 30 (blue) and L = 50 (black). 
Averaged over 5000, 500 and 50 runs respectively. 
Data used for table 3. 



